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Pair fluctuation theory has been used to study the crossover from the weak coupling BCS theory 
to the strong couphng Bose Einstein Condensation. The effect of fluctuations has been studied over 
the whole crossover regime. It has been shown that the pair fluctuations are enhanced considerably 
in both two and three dimensions and hence mean fleld theory is inadequate to study the physical 
properties in this regime. A self consistent scheme for calculating the pair susceptibility is given. 



The study of crossover from BCS to Bose Einstein Condensation has been a subject of intense study in recent 
times. The motivation for this study has been from the point of view of trying to model High Tc superconductors. 
These materials share many properties with BCS superconductors, however, the coherence length ^ in these materials 
is unusually small [n|. The short coherence length feature is also apparent in the phenomenon of Bose Einstein 
Condensation (BECj where the coherence length ^ is expected to be of the order of an atomic spacing. Moreover, 
in the BEC process the temperature at which pairing occurs is different from the pair condensation temperature. 
The phenomenon of BEC was was known as Schafroth condensation |^ in the pre BCS days of superconductivity. 
Recently this has lead to a formulation of theories based on occurance of pseudo gaps in the normal state of these 
materials. [^|-^The crossover from BCS to BEC has been studied as a case of weak to strong attractive coupling 
theory in details during recent years within the mean field theory or RPA [0-0. These results agree with the 
expected behavior at the two extreme limits (i.e. the weak coupling BCS limit and the strong coupling BEC limit) 
and predict a smooth crossover from one limit to the other. However, the intermediate regime, to which the High 
Tc superconductors are believed to belong within this theory, has not been well understood. The theory is based on 
the existence of two parameters which can be varied - the electronic density n and the coupling strength g. Coupled 
self consistent equations for these two parameters leads to the calculation of Tc, the transition temperature and the 
chemical potential /i. For a constant density, if the coupling strength is varied, there is smooth crossover from the 
weak coupling BCS superconductivity to the strong coupling Bose Einstein condensation. In principle, however, both 
the parameters can be varied. If this is done such that Tc remains constant and fi is varied there is a crossover with 
respect to fj,, where for large positive fj, {p ep, the Fermi temperature) the results for the BCS superconductors are 
recovered and for large negative fi the BEC results are recovered. Thus, the intermediate regime can be interpreted 
in this regime as the region close to /i = 0. In what follows, we shall take this point of view in order to study the 
crossover. 

In this work, we wish to study the crossover regime from the point of view of pair fluctuations. The main aim of this 
work is to ask the question - how important are fluctuation contributions to this phenomenon of crossover, especially 
in the region close to the crossover ? As has already been mentioned before, most of the studies in this subject have 
taken recourse to mean field theory. If fluctuation contributions do become important, then, mean field theory is 
bound to breakdown. We therefore set up an equation for the pair susceptibility within the mean fluctuation field 
approximation (MFFA) where the effect of fluctuations are included in a self consistent manner. This self consistent 
equation is then solved to obtain the inverse pair susceptibility over the whole regime of crossover from the BCS limit 
(where ^ > T) to the BEC hmit (where > T). 

The starting point for the calculation is the BCS Hamiltonian in the momentum space with an s-wave attractive 
potential which is given by, 

H — — epCj, g.Cp,a — —g Cp^g ,j.C^_p _„C-p'^-crCp'+q^a, (1) 

P;<7 p,p',q,(7 

where ep = /2ra is the free electron kinetic energy, and g is the interaction strength. The mean field results for this 
Hamiltonian are well known for both the BCS limit and the BEC limit. 

The functional integration technique is a convenient method for treating the fluctuation problem microscopically 
0. For the case of superconductors the technique was developed by Langer Rice jl^ and others. We follow 
the procedure adopted by Rice in what follows. Using the functional integral technique the partition function can be 
expressed as, 

Z/Z, = y n exp (2) 
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where is the partition function of non-interacting system and Xq^m is the (bosonic) auxihary field with momentum q 
and frequency uim- The pairing correlation is the only one which has been considered here. A diagrammatic expansion 
of F[a;g,m] in powers of Xq,m is possible and leads to the Ginzburg Landau functional. The corresponding expression 
for the free energy functional is given by, 

-F'' ''[s^g.in] = ■^^^(1 + '^)^q,m^q,m + "^3 -^qi ,''(32 ,93 '^92 ■^93 -^ijl +«2 -<23 ■ (3) 

q,m 9,^ 

These terms are shown in Fig. (||). 
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FIG. 1. Lowest order diagrams for Free Energy 

The Eq. |^ represents the free energy of interacting electrons near a superconductive instability in terms of interacting 
pair fluctuation field. The first term is the "free" or the RPA term (Fig. |l|a), while the second term represents the 
interaction between pair fields (Fig. ^3). Parameters of this model, e.g. transition temperature, collective mode 
dispersion, fluctuation spectrum, fluctuation coupling vertices are determined by properties of the underlying fermion 
system. Expanding K^^^^ and K^^') in powers of q and uj„i and retaining the lowest powers we get. 



PF[Xq,m] = [a + Cq'^ - lUJyn] \ Xq.m P +b \ Xq^rn I'' 

The coefficients a, b, c and d are given by. 
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where vi — {21 + 1)tt/(3. The value of comes out to be 2 in 2D and 4/3 in 3D. At T = Tc, a = 0, which leads to the 
Thouless criterion, (1 + qKq = 0), 
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tanh[/3,(e-/ie)/2] 
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The integral on the right hand side is UV divergent. In the conventional weak coupling theory this is taken care of by 
putting an upper cutoff at the Debye energy. However, since we are interested in the crossover from weak to strong 
coupling behavior this restriction on the energy integration has to be removed. The singular behavior is taken care 
of by replacing the bare g by the low energy limit of the two body T matrix. 
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where 1/t turns out to be m/ Alias for d = 3 and (m/47r) ln(i?b/2) for d = 2. Here as is the s-wave scattering length 
and Ei, = 1/ma^. 

The number density of electrons n and the coupling strength g determine the transition temperature and also 
specify whether one is in Bose or the BCS regime of pairing. These are related to the transition temperature Tc and 
the chemical potential /i through Eqs. ^, |^ and the equation for the number density, viz. 

.^^(1-tanh^i^), (8) 

k 

which are to be solved self consistently. One can consider [n and g) or, (Tc and ^) combination as the parameters in 
the problem. We find convenient to consider the latter pair, particularly when dealing with the crossover region. The 
region around /i ^ marks the crossover. We write a etc. in terms of Tc and , i-c, 

a='-(r ,,,(,)tanh[/^.(.~M.)/2] _ r tanh[/^(. ,)/2]N 

2 (e - Mc) 7o (e - / 

The expressions for 6, c and d are nonsingular at Tc ■ 

The pair susceptibility can be calculated in this scheme as, 

X{q,UJm) =<\ Xq,m P>= j dXq^rndx*^^^ \ Xq^rn P CXp -/3i^[a;q,,„]/ j dXq^„idx*„^exp- PF[Xq^rn]- (10) 

We now write down a self consistent equation for the inverse of pair susceptibility within the mean fluctuation field 
approximation. This amounts to writing the quartic term, in the expansion of the free energy functional in terms of 
the pair field, as a quadratic form (i.e. | Xo,m l^^^^l 2^9, m P>| Xq^m p) and taking average over the effective Gaussian 
distribution. The technique is well known [|17| ; we, however, follow the earlier work on spin fluctuations 

111- The 

pair susceptibility is given by, 
where, 

a{T,fi)^a + b ^ x{q',^^ni')- (12) 




FIG. 2. Mean fluctuation fleld diagram 

To calculate the susceptibility one needs values of Ginzburg-Landau coefficients, a, b, c and d. On performing the 
summation over the Matsubara frequencies, one gets in 3 dimensions. 



^ = o / de p3n(e)$(e), c = — / de P3d{^)(- $(e), and d = - / de psD 
° Jo Jo o Jo 



(13) 



Here, P3D(e) — (2m)'^/^Y^/47r^ is the density of states in 3 dimensions and 

$(e) = {2tanh(/3(e - /i)/2) - /3(e - M)sech2(/3c(e - ^^)/2)}/{e - (14) 

It is not possible to evaluate the energy integral analytically throughout the parameter space. Only in the extreme 
limits the results are available. In the BCS limit Pfi ^ 1, a = Inije^ /Se^ (3ijl)p3d{p)) — ln(7re^/8e'^/3/Xc)P3D(Mc), 
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b = (7C(3)/87r2)p3£,(/x)/52 c = {7C{3)/12^^'^)p3D(p)^J■P'^, and d = ln(8eT/3/i/7r)p3r,(/^)/V ■ In the Bose Einstein 
regime {-(3fi > 1) the corresponding value are, a = {tt/2){p3d{\ M I) - PsdH IJ-c |)), b = P3d{\ p |)(7r/32 | p p), 
c = /33d(| M |)(7r/16 I pi I), and d = p3i5(| ^ |)(7r/8 I I)- 

Similarly, in two dimensions, the summation over the Matsubara frequencies gives, 

, ^ / ^ 1/^7C(3)^2 ^ , , , , tanh(/3e/2) \ ^ , , , tanhifi p / 2) 
b^^J^ dep2D{eMe), c ^ - l^^P^p9{p)+ \ p \ ^ d^p^oie) 1 ' ^'^^^ d-p2D(e) ^f^' 

(15) 

where p2D(e) = m/27r. Again, in the BCS limit, a = P2d\i^{I3cPc/ fip), b = (7C(3)/87^2)p2D/3^ c = (7C(3)/87r2)p2DAi/3^, 
and d — P2_d/4/^; and in the Bose Einstein limit, a — (l/2)/92Dln(/^)/^c, b ~ P2L>/8/i^, c = P2_d/8 | | and d = P2n/4 | 

Since the values for the GL coefficients are expressible analytically only in the two limiting cases we use the full 
expression for these coefficients to calculate the susceptibility and take recourse to the numerical evaluation of the 
integrals. 

Before proceeding further we want to make some remarks regarding these coefficients. In the BCS limit the 
temperature dependence of c (which is related to and the coefficient of the fourth order term are similar. (In fact 
b and c are related in both the limits). The reason for this can be traced to the diagrammatic representation of these 
terms. The quadratic term has two internal fermion lines and to get the coefficient c we differentiate K^^^ twice with 
respect to momentum. This makes it equivalent to four fermion lines which is similar to four fermion lines in the fourth 
order term, particularly when vanishing momentum limit is taken. Moreover, b and c are singular as T ~ Tc —> 0. 
This may have interesting consequences when one considers the possibility of quantum phase transition in this limit. 
The second remark is regarding the BE limit. In this limit the coefficients, in both two and three dimensions, are 
independent of temperature. The GL expansion look similar to one for system near an antiferromagnetic instability 
[ p^ . Thus in this limit, one can expect that in the limit of Tc ^ (in case it is possible) a non- Fermi liquid like 
behavior in some transport and thermodynamic properties may be seen and the temperature dependence of these 
quantities will be similar to that in antiferromagnets. 

To calculate the pair susceptibility the equation for a{T,p) must be solved self-consistently. For this we first note 
that, 

1 >r- . ^-i r , ^ 

P^ ^ ^^ ^ ^ ^J-oo i^iT) + cq^r + dP^^ (e^- - 1) 

= ^[lny-(2y)-i-V(2/)], (16) 

where i/'(2/) is the digamma function and y is given by, (a + cq^) /2T:dT . A good interpolation for In y — {2y)~^ — 'ipiu) 
which is valid for large as well as small y is l/(2j/+12?/^). Thus, the self-consistent equation for a{T,p) can be written 
as, 

aiT,p)=a + -J^ dq^^^^, (17) 

where D denotes the dimension of the system. Our interest is in how a(r, p) which has been obtained within the 
fluctuation theory compares with the mean field value a. We calculate a(T, p) as a function of p at temperatures just 
above Tc. 

The results of self consistent calculation of Eq. ( pT[ ) with the numerical evaluation of the coefficients 6, c and d 
in three and two dimensions is plotted along with the corresponding mean field value a in Fig (^). The parameters 
are taken as (3c = 100 and (3 = 90. The behavior of the curves is not very different in two and three dimensions 
qualitatively. In the calculation of the coefficient a, we have assumed that p at temperature just above Tc is the same 
as its value at Tc, i.e. pc- This is reasonable for temperatures close to Tc- This is why in the extreme negative p 
regime (the BEG regime) a is identically equal to zero since in this regime a is independent of temperature. From 
the figure two things stand out very clearly - the first is that in the extreme limits (both the BGS and the BEG), 
the effect of fluctuations is negligible. Secondly, as one moves away from the two extreme limits by reducing | |, 
the effect of fluctuations increases leading to a dominant contribution to the pair susceptibility as /i — > 0. In the the 
BEG regime as one moves towards the crossover regime the value of a also increases but the fluctuation contribution 
is much larger. It is also apparent that, the fluctuation contribution in two dimensions is an order of magnitude 
higher than in three dimensions for the same set of parameters. The curves are shown for relatively small transition 
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temperature. Preliminary results indicate that as the parameter /3c is reduced from the value 13c = 100 towards 
(3c = 20 corresponding to an increase in the transition temperature, the crossover region as well as the fluctuation 
contribution both increase. 
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FIG. 3. Inverse susceptibility within the mean fluctuation field approximation {a(T,n)) and within the mean field approxi- 
mation (a(r, /x))plotted as a function of /x in three dimensions (a) and in two dimensions (b). Note that the scale for the plot 
in two dimensions (b) is one order of magnitude larger than that in three dimensions(a). 



In conclusion, we have studied the crossover from BCS to BEC from the point of view of pair fluctuations. We find 
that the effects of fluctuations are large in the crossover regime, where /i is close to zero. This shows that a mean field 
calculation will be inadequate in this regime. The enhanced fluctuation effects shows that the behavior of physical 
properties can be qualitatively different in this regime. If a parallel is drawn with the results of Spin Fluctuation 
theory, then it can be predicted that the physical properties of the system may show a non-Fermi liquid behavior if 
parameters are such that one is close to Tc with Tc — > 0. The non-Fermi liquid behavior of some high Tc materials 
can be modeled in this scheme, our discussion has been confined to isotropic two- or three- dimensional systems, 
while the anisotropy seem to play an important role in the High Tc materials. This is especially true for the layered 
High Tc compounds. However, in the layered materials, the coherence length along the a—b plane is also of the order 
of that expected for the bulk material. Hence, the normal state properties of the a — b plane, which are different 
from the properties along the c axis, should in fact follow our two dimensional result. The other important point is 
that the form of the pairing potential which we have chosen here is s-wave. It is now quite well agreed upon that 
the pairing potential in most of the High Tc materials, especially the layered materials, is d-wave. In our opinion the 
enhancement of fluctuations will occur in this case also along with the reflection of the inherent anisotropy brought 
in by the d-wave nature of the pair wave function. Work in this direction is in progress. 
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